Cosmological viability of f{R) - gravity as an ideal fluid and its compatibility with a 

matter dominated phase 
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We show that /(_R)-gravity can, in general, give rise to cosmological viable models compatible 
with a matter-dominated epoch evolving into a late accelerated phase. We discuss the various rep- 
resentations of /(i?)-gravity as an ideal fluid or a scalar-tensor gravity theory, taking into account 
conformal transformations. We point out that mathematical equivalence does not correspond, in 
several cases, to the physical equivalence of Jordan frame and Einstein frame. Finally, we show that 
wide classes of ,f{R) gravity models, including matter and accelerated phases, can be phenomenolog- 
ically reconstructed by means of observational data. In principle, any popular quintessence models 
could be "refrained" as an /(_R)-gravity model. 
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1. Up to ten years ago, the paradigm of cosmology has 
been the Cosmological Standard Model inspired by Ein- 
stein's General Relativity (GR) , the Big Bang model and 
particle physics. Recently, several observational surveys 
have provided high quality data giving rise to the so 
called precision cosmology by which the nowadays picture 
of the universe results radically different from the stan- 
dard Einstein-Friedman model. Specifically, the Hub- 
ble diagram of Type la Supernovae 1], anisotropics in 
the cosmic microwave background radiation 2], matter 
power-spectrum measured in optical surveys of large scale 
structures 3J have convincingly given the picture of a 
spatially flat universe in a phase of accelerated expansion 
(for a review see |j|). The new theoretical scenario seems 
consistent only if huge amounts of dark componentSjUSu- 
ally treated as fluids, are introduced into the game |^|^. 
In particular, in order to drive the observed cosmic accel- 
eration and to fill the gap between the matter-energy con- 
tent and the critical density of a spatially flat universe, 
a negative pressure fluid, not clustered at small scales, is 
needed. The dramatic result is that this mysterious and 
unexpected component, referred to as Dark Energy, adds 
to the already supposed presence of the Dark Matter nec- 
essary to fit the astrophysical data at smaller clustered 
scales. In the simplest scenario (Concordance A-CDM 
model) Q, dark energy is interpreted as the cosmologi- 
cal constant, contributing for 70% to the whole energy 
budget, while the remaining 30% is constituted for 4% 
of baryons and for 25% of cold dark matter. Such a kind 
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of dark matter is typically addressed to exotic particles 
such as, e.g., WIMPs or axions, nevertheless, up to now, 
no definitive proof of their existence has been achieved 
■8|. This model, although satisfactory from an observa- 
tional point of view, is anyway theoretically disfavored 
since there is a huge (~120 orders of magnitude) differ- 
ence between the predicted and observed values of A. 
In order to overcome the cosmological constant problem, 
many alternative models have been proposed. Unfortu- 
nately, none of these attempts is fully satisfactory, since 
they mostly lead to 'ad hoc' cosmologies, not well the- 
oretically and/or observationally founded (for a review 
seeHIIg). 

A different approach to face this challenge can be pur- 
sued reversing the problem and considering extended 
modified gravity compatible with the observational re- 
sults. This scheme possesses the relevant feature that 
experimental data turn out to be naturally interpreted 
without the need of additional components/fluids. Of 
course, modifying the gravity from the Einstein GR im- 
plies several theoretical issues. In particular, one has 
to be able to match GR prescriptions, at least at the 
post-parameterized- Newtonian limit 'll', and recover the 
positive and well-established results of standard cosmol- 
ogy. The attempt of extending GR is indeed theoret- 
ically and physically appealing. In fact, rather than 
formulating a different approach radically changing the 
rules on the standard one, one can straightforwardly ex- 
tend the properties of Einstein gravity. In such a way, 
the theoretical foundations remain the same, but the 
search for a different theory of gravity can be essentially 
pursued by inferring physical features directly from the 
data. Actually, among Extended Theories of Gravity, 
/(i?)- gravity [lOl Il2. Il3l | represents a viable alternative 
to dark energy and naturally gives rise to accelerating 
sing ularity- free solutions in early and late cosmic epochs. 
HJi HM HM H^ nil ' Furthermore, it has to be remarked 



that higher-order terms in the gravity Lagrangian are re- 
quired in theories attempting to formulate a satisfactory 
quantum field theory on curved space-times [l2| and can 
descend from the low energy limit of String/M-theory 
|20] . These theories fit quite well the data on the SNela 
Hubble diagram |22] and can also give interesting theoret- 
ical predictions with respect to the CMBR observational 
results 113. 

f{R) Lagrangians in the gravitational action, con- 
ceived as analytic functions of the Ricci scalar, deserve a 
special attention. These actions are the simplest among 
those which can be constructed by the curvature tensor 
invariants and their covariant derivatives ^; they allow 
to determine stable models with respect to small pertur- 
bations around maximally symmetric solutions and to 
avoid ghost-like behaviors in effective field theory. Be- 
side cosmology, such models could have interesting con- 
sequences at galactic and Solar System scales. In fact, al- 
though still controversial, the PPN limit of several f{R) 
theories of gravity could be reconciled with the gravita- 
tional experiments, furnishing a theoretical background 
for some measured deviation from GR ,23_^ 2^. How- 
ever, these results need a deeper investigation and such 
a purpose will be pursued together some gravity ded- 
icated experiments |25j . However it is worth stressing 
that f{R) theories, in particular their simplest power law 
form f{R) = /o-R", represent a viable model even from 
the astrophysical point of view. It is in fact possible to 
demonstrate that they can furnish a natural mechanism 
to interpret the astrophysical data at galactic scales. Ac- 
tually, thanks to the corrections to the Newton potential 
which they give rise in the low energy limit |42ll43l| . one 
could achieve a geometrical explanation to the missing 
matter problem (as suggested also by other alternative 
approaches as MOND [23 , which recently have been re- 
lated to higher order gravity itself |27j . 

An important issue regarding fourth order/modified 
gravity is represented by its peculiarity to be recast in a 
(non) - minimally coupled scalar-tensor theory by means 
of a conformal transformation, namely, from the Jordan 
frame to the Einstein one. In this way, in principle, one 
can pass from higher order gravity to scalar-tensor grav- 
ity and vice-versa. From a formal point of view, this 
procedure is absolutely plausible while its physical rele- 
vance is still controversial, [2^ |23, 113 and, up to now, 
there is no privileged point of view regarding the physi- 
cal relevance of conformal transformations which remains 
still debatable ISOJ . 

Following standard approaches, one cannot conclude 
univocally about the physical relevance of the results. In 
general, it is customary to pursue conformal transforma- 
tions with the aim of studying a physical problem de- 
veloped in the Jordan frame and then mapping it in the 
Einstein frame where the equations turn out to be math- 
ematically simpler. Therefore, Einstein frame equations 
are solved and the solutions are considered as physically 
significant [3l|, Hj, 113 ■ Despite several results achieved 
along this line, in order to get the physical meaning of 



the problem, a back mapping of the results obtained 
in the Einstein frame to the original Jordan one (see 
i.e. |34L II3) has to be in order. In general, the back- 
transformation to the Jordan frame does not assure to 
get well endowed physical solutions. It has been even 
widely demonstrated that passing from one frame to the 
other can completely change the physical meaning as well 
as the stability of the solutions one is dealing with 3g| . 
For example a scalar field quintessence model can roll a 
V{(f>) potential motivated by particle physics in the Ein- 
stein frame, while its back-transformation to the Jordan 
frame corresponds to a potential V{(j)) with no physical 
interest. Furthermore, in the same sense, if a solution 
is accelerating in one frame, its conformally transformed 
counterpart in the other frame is not necessarily acceler- 
ating too 33] , or equivalently, if one obtains a phantom 
model in the Einstein scheme, its mapping in the Jordan 
frame does not give, in general, a phantom-like evolution 
[37l |. For example, if one takes into account a conformally 
coupled scalar field with the potential V{(l>) = A^"* in the 
Einstein frame '33] , the conformal anti-transformation of 
such a potential gives : 
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which is difficult to be physically interpreted. Such a 
problem is particularly relevant if standard matter is con- 
sidered in the field equations due to non-minimal cou- 
plings arising thanks to the conformal transformations. 
As an example, there is a recent result where a solution 
obtained from /(i?)-gravity, transformed in the Einstein 
frame and connecting the radiation epoch to the acceler- 
ating attractor, a ~ t^'^ , is considered non-physically rel- 
evant since its back transformation in the Jordan frame 
gives a{t) — t^/'^ which does not accord with the ex- 
pected dust matter dominated regime [l^- The inter- 
pretation of the result is difficult and controversial since 
assuming a priori the Einstein frame as " physical" , we 
should rule out /(i?)-theories. This means that, despite 
several positive results achieved from extended modified 
theories of gravity, without some hints coming from the 
original framework in which the theory is conceived there 
is no way for a clear physical interpretation of what is go- 
ing on. In other words, is it correct to obtain a solution 
in a frame and then interpret it in another frame? 

In this paper, we discuss how extended gravity (in par- 
ticular /(i?)-gravity) can be dealt under the standard of 
an ideal fluid description or transformed to a (mathe- 
matically) equivalent scalar-tensor theory. We show that 
the three descriptions could not be physically equivalent. 
For example, the phantom phase in one of the three ap- 
proaches can corresponds to non-phantom phase into an- 
other one. Furthermore, it is possible to show that a dust- 
matter regime can be obtained in a /(i?)-gravity sce- 
nario, and, in particular, its evolution can be connected 
to a late time stage of accelerating expansion, without 
considering conformal transformations. In Appendix it 



is demonstrated how any realistic cosmology may be re- 
alized via the reconstruction of modified gravity. These 
facts point out that the mathematical equivalence among 
Jordan, Einstein and "fluid" descriptions does not neces- 
sarily implies physical equivalence and solutions should 
be carefully studied into the frames in which they are 
obtained. 

2. In this section, we show that modified gravity can 
be transformed also to GR with standard ideal fluid, so 
actually two Einstein frames exist! 

Let us start from the action of the /(i?)-gravity: 
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Here f{R) is a function of the scalar curvature R. For 
example. 
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is a generalization of the model in [ig . In Q , as long as 
n > —1, n is not restricted to be an integer but can be 
an arbitrary real number. As a more generalized model, 
we may consider [l7j 
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Here we assume to > 1. When the curvature is small, 
the behavior of the model Q) is not changed much from 
that of Q but the model may describe the inflation in 
the early universe. In the following, we assume ci^2 > 0. 
For the model lj21l, Starobinsky has proposed a condi- 
tion 
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which could be necessary for the stability of FRW solu- 
tions and for the existence of the positive mass squared 
for the "scalaron". For the model Q, the condition gives 



n{n + l)ci 
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which could be satisfied if — 1 < n < when ci > and 
R> 0. And for the model l@J, we find 



n(n + l)ci 



i?"+2 
Then if ci, C2, i? > 0, wc find 

n(n + l)ci 
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The condition is satisfied again if —1 < n < 0. When 
n > (we do not consider the case n = since it corre- 
sponds to the cosmological constant), the condition © 
is satisfied if C2 is large enough. Since the condition © 



depends on the curvature, even if the condition could be 
satisfied in the present universe, when the curvature be- 
comes smaller, the condition may be violated. This 
may indicate that even more terms should be considered 
in modified gravity action in the future universe. Cur- 
rently, such terms may be next-to-leading order but they 
may become essential at extremely small curvature. 

Technically, the action Q can be rewritten in the form 
of the scalar-tensor theory. By introducing the auxiliary 
fields, A and B, one can rewrite the action j^Jl as follows: 



S^ d^x^^{B {R-A) + f{A)} . (9) 

One is able to eliminate B, and to obtain 

S^ J d^x^^ifiA) iR-A) + f{A)} , (10) 

and by using the conformal transformation g^^ — + e'^g^i, 
(cr = — In (2k^ f (A)^) , the action pO|l is rewritten as the 
Einstein-frame action: 
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and G is defined by G (2k^ f (A)^ = A, that is the inverse 
function of 2k^ f . In case of (O, we find 



U{a) = {n+l){2K^a) 



-1 



(13) 



In principle, the inverse process is possible, that is, by 
starting from the action of the Einstein gravity coupled 
with scalar field, one can obtain the f{R) action. By 
writing 
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we can start from the following action: 
S= I d^x 



g^^R-^d^vd'^^-Vi^) 



(15) 



Then one can use the back conformal transformation 
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and make the kinetic term of (p vanish. Hence, one ob- 
tains 



S = / d^xy^ ■ 
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-R 



^±2Kip 
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Since Lp becomes the auxiliary field, one can delete Lp by 
using as an equation of motion: 



R = e±«^Vl Un^V{^) ± 2k^-V'{^)\ 



(18) 



which can be solved with respect to ip as ip — 'p{R). The 
action (|17|l is re-incarnated in the form of /(i?)-gravity, 
that is: 



f{R) 
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Therefore /(i?)-gravity theory is mathematically equiva- 
lent to the scalar-tensor theory which does not yet mean 
their physical equivalence |37|. 

Let W be the inverse function of V{(p): W{V{(p)) — (p. 
Then if we can solve the following differential equation^ 
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(20) 
with respect to /i((/'), we can further rewrite the action 
(Uni as 



S = I d^x^—g { ^R - ^uj{cb)d^^d'^cf, - V{cf>) 



(21) 



Here 
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and 
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For the action (|21|l with (|22|l , if we can neglect contribu- 
tions from other matter fields, the following solution of 
the FRW equations can be found |4(TI |: 



t , H ^ h{t) 



(24) 



Here p and p are the scalar energy density and the pres- 
sure defined as: 



1 



2 + V{cj>) , p=-u; 



y(0) . (26) 



We should note that there is almost one to one (up to 
integration constants) correspondence between V{lp) in 
(|22|l and /i((/)), and therefore the corresponding EoS ideal 
fluid (|25|l . Then for arbitrary /(_R)-gravity, in principle, 
one can find the corresponding ideal fluid EoS description 
(what maybe considered as yet another Einstein frame). 
As an example, the small curvature case for (jSJ may 
be considered. Hence, 
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In (|28|l . the two signs ± are independent with each other. 
Then W (HOJ has the following form: 
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Here x is an a-dimensional variable, which is used to 
express the function W . Then a non-trivial solution for 
(Pn|l is given by 



/i((/)) = 3 
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Using (|25(l . we find an EoS with constant EoS parameter 
w: 



n,p. -=9^^) -1- (31) 



For instance, if 
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Furthermore one finds the following, effective equation of 

state (EoS): 



we find w = 0, which corresponds to dust. On the other 
hand, if 
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^ A trivial solution of l2Ui is h{<l>) = ho with a constant hg. The 
solution corresponds to de Sitter universe. We should note that 
/(-R)-gravity, in general, admits the de Sitter solution. The h = 
ho solution could correspond to the de Sitter solution in f{R)- 
gravity. 



it follows w — 1/3, which corresponds to radiation. 
Hence, extended modified gravity may be presented in 
mathematically equivalent form as GR with ideal fluid! 

Note that as we have used the scale transformation, 
the cosmology appears in a different way in the frame 
corresponding to (O from the frame corresponding to 



(|15|l . In the frame corresponding to H15() . which we caU 
the Einstein frame, by using (|24|) and H3U|) . we find 



H = 3 



or a oc i 



■^ln + 2 j 



and by using (O, (|2S1), and (fT?)|l . 
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Here io is a constant of the integration in H23I) . The 
time coordinate ij and the scale factor aj in the frame 
corresponding to ©, which we call as the Jordan frame, 
is related with t and a in the Einstein frame as 



dti^c'^/^dtoit^i'^Tdt 
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Hence, one gets 
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The latter Ea. (|37|l reproduces the standard result. We 
should note that the cosmology is really changed in the 
two frames. In the Einstein frame H34I) . we find always 
H < 0, which corresponds to non-phantom phase. On 
the other hand, in the Jordan frame, we find 



_ aj _ {n + l)(2n + l) 1 



aj 



t 
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Then if — 1 < n < —1/2, it follows H > 0, which corre- 
sponds to phantom phase. We should note that in the 
phantom regime —1 < n < —1/2, Starobinsky's condi- 
tion © is satisfied. Hence, we demonstrated that math- 
ematical equivalence does not always means the physi- 
cal equivalence (the same point of view was recently ex- 
pressed in [41|). In conclusion, we can state that the 
existence of so many different frames indicates that only 
one, the original frame, can be considered the physical 
one. 

3. As we have seen, by performing a conformal trans- 
formation on a modified gravity Lagrangian f{R), it is 
possible to achieve, in the Einstein frame, dust matter 
behaviours which are compatible with observational pre- 
scriptions. In addition, by exploiting the analogy be- 
tween the two frames and between modified gravity and 
scalar-tensor gravity one can realize that physical re- 
sults, in the two conformally related frames, could be 
completely different. In other words one can pass from 
a non - phantom phase behaviour (Einstein frame) to a 
phantom regime (Jordan frame). Now, we can suppose 
to change completely the point of view. In fact, we can 
rely directly with the Jordan frame and we can verify 



if a dust matter regime is intrinsically compatible with 
modified gravity. 

As a first example, one can cite the exact solution pro- 
vided in 21], which has been deduced working solely in 
the Jordan Frame (FRW universe). In particular, one is 
able to find a power law regime for the scale factor whose 
rate is connected with the power n of the Lagrangian 
f{R) — foR"- In other words, one has a{t) = agP with 
^ _ -2n+Zn-\ ^ gy(,j^ g^Q exact solution is found out 
when only baryonic matter is considered |42l l43l | . It is 
evident that such a solution allows to obtain an ordinary 
matter behaviour (a = 2/3) for given values of the pa- 
rameter n (i.e. n ~ —0.13, n ~ 1.29). Such solutions 
are nevertheless stable and no transition to acceleration 
phase then occurs. In general, it is possible to show that 
solutions of the type 



ao(t-to)='<^+™' 



(39) 



where w is the barotropic index of standard perfect fiuid, 
arises as a transient phase, and this phase evolves into an 
accelerated solution representing an attractor for the sys- 
tem [ij . In any case, a single solution exactly matching, 
in sequence, radiation, matter and accelerated phases is 
unrealistic to be found out in the framework of simple 
/(i?)-power law theories. The discussion can be further 
extended. Modified gravity can span a wide range of an- 
alytic functions of the Ricci scalar where /(i?) = /o-R" 
only represents the simplest choice. In general, one can 
reverse the perspective and try to derive the form of grav- 
ity Lagrangian directly from the data or mimicking other 
cosmological models. Such an approach has been de- 
veloped in 45] , and allows to recover modified gravity 
Lagrangians by the Hubble flow dynamics H{z): in par- 
ticular, it is possible to show that wide classes of Dark 
Energy models worked out in the Einstein frame can be 
consistently reproduced by /(i?)-gravity as quintessence 
models with exponential potential [461] . Clearly the ap- 
proach works also for the case of coupled quintessence 
scalar field. In other words, the dynamics of H{z), con- 
sidered in the Jordan frame, is reconstructed by observa- 
tional data considered in the Einstein frame then assum- 
ing one of the two frames as the "physical frame" could 
be misleading. Here we further develop this approach 
with the aim to show the viability of f{R) gravity to re- 
cover a matter-dominated phase capable of evolving in a 
late accelerating phase. 

From a formal point of view, the reconstruction of the 
gravity Lagrangian from data is based on the relation 
which expresses the Ricci scalar in terms of the Hubble 
parameter : 



R = -qIh + 2H^ + —J . (40) 

Now, the modified gravity field equation reads [isl |21| : 



Gap - Rap - -^RgaP - T^p"""" + T^ (41) 



where a stress - energy tensor summing up all the higher 
order contributions is defined. It plays the role of an 
effective curvature fluid : 
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+ nR)''''''{ga^.9pu-gc.pg^u)} (42) 

while the matter term enters Ea. (|41(l through the modi- 

fled stress -energy tensor T^™' = f^^^/f{R) with f^"^ 
the standard minimally coupled matter stress - energy 
tensor. Starting from this general scheme, one can re- 
construct the form of f{R) from the Hubble parameter 
as a function of the redshift z exploiting the relation H4l)|l 
after this expression has been rewritten in term of the 
redshift itself. A key role in this discussion is played by 
the conservation equation for the curvature and the mat- 
ter fluids which, in the case of dust matter, (i.e. p„i = 0) 
gives : 



Pcurv ~r <J-i^ t^-L n^ ^curv ) Pc 



-J77^{Pm+3Hp,rr) 
dfiR) ,,„. 



In particular, one may assume that the matter energy 
density is conserved : 

Pm = nMPcr^ta'^ = 3H^nM{l + zf (44) 

with z = l/a — 1 the redshift (having set a(io) — 1), ^m 
the matter density parameter (here and hereon quantities 
labelled with the subscript refers to present day {z = 0) 
values). Ea. l44|) inserted into Ea. (|43|) . allows to write a 
conservation equation for the effective curvature fluid : 
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[riRW 

Actually, since the continuity equation and the fleld equa- 
tions are not independent |43 , one can reduce to the fol- 
lowing single equation 
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+R 



L-^{3H^nM{l + zf + Rr{R) + 
Rf'iR) - Hf"{R)] } , 
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where all quantities can be expressed in term of redshift 
by means of the relation ^ = —{l + z)H-j^. In particular, 
for a flat FRW metric, one has : 
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Now, we have all the ingredients to reconstruct the shape 
of f{R) by data or, in general, by the deflnition of a suit- 
able H{z) viable with respect to observational results. In 
particular, we can show that a standard matter regime 
(necessary to cluster large scale structure) can arise, in 
this scheme, before the accelerating phase arises as in the 
so called quiessence model. A quiessence model is based 
on an ordinary matter fluid plus a cosmological compo- 
nent whose equation of state w is constant but can scat- 
ter from w = — 1. This approach represents the easiest 
generalization of the cosmological constant model, and 
it has been successfully tested against the SNela Hubble 
diagram and the CMBR anisotropy spectrum so that it 
allows to severely constraint the barotropic index w |47l| . 
It is worth noticing that these constraints extend into 
the region w < —1 therefore models (phantom models) 
violating the weak energy condition are allowed. From 
the cosmological dynamics viewpoint, such a model, by 
definition, has to display an evolutionary rate of expan- 
sion which moves from the standard matter regime to 
the accelerated behaviour in relation to the value of w. 
In particular, this quantity parameterizes the transition 
point to the accelerated epoch. Actually, if it is pos- 
sible to find out a /(it!)-gravity model compatible with 
the evolution of the Hubble parameter of the quiessence 
model, this result suggests that modified gravity is com- 
patible with a phase of standard matter domination. To 
be precise, let us consider the Hubble flow deflned by this 
model : 



H{z) - Ho^nMil + z^ + nx{l + z)3(i+-) (51) 

with fix — (I — Qm) and w the constant parameter defln- 
ing the dark energy barotropic index. This deflnition of 
the Hubble parameter implies^ : 

R = -ml klMil + zf + nx{l~ 3u;)(l -f z)3(i+-) . 

(52) 
The ansatz in Eq.|(5lJ allows to obtain from Ea. (|46|) a 
differential relation for f{R(z)) which can be solved nu- 
merically by choosing suitable boundary conditions. In 
particular we choose : 



dR 

dz 



(53) 



/■<«'H'f)'| 



(48) 



^ Note that R is always negative as a consequence of the signature 
{-f-, — , — , — } adopted. If we had used the opposite signature, 
Ea. l52i is the same, but with an overall positive sign. 









f{z = 0) = /(i?o) = 6i/2(i - ^m) + Ro ■ 



(54) 



(55) 



A comment is in order here. We have derived the present 
day values of df /dz and d^ f /dz^ by imposing the consis- 
tency of the reconstructed f{R) theory with local Solar 
System tests. One could wonder whether tests on local 
scales could be used to set the boundary conditions for a 
cosmological problem. It is easy to see that this is indeed 
meaningful. Actually, the isotropy and homogeneity of 
the universe ensure that the present day value of a what- 
ever cosmological quantity does not depend on where the 
observer is. As a consequence, hypothetical observers liv- 
ing in the Andromeda galaxy and testing gravity in his 
planetary system should get the same results. As such, 
the present day values of df /dz and d'^f/dz'^ adopted 
by these hypothetical observers are the same as those 
we have used based on our Solar System experiments. 
Therefore, there is no systematic error induced by our 
method of setting the boundary conditions. 

Once one has obtained the numerical solution for f{z), 
inverting again numerically Ea. (|52|l . we may obtain z = 
z{R) and finally get f{R) for several values of w. It 
turns out that f{R) is the same for different models for 
low values of R and hence of z. This is a consequence 
of the well known degeneracy among different quiessence 
models at low z that, in the standard analysis, leads to 
large uncertainties on w. This is reflected in the shape of 
the reconstructed f{R) that is almost w - independent in 
this redshift range. An analytic representation of the re- 
constructed fourth order gravity model, can be obtained 
considering that the following empirical function 



ln( 



-/) - h [In {-R)t [1 + In {-R)t + h (56) 



approximates very well the numerical solution, provided 
that the parameters (^i, Z2, ^3, ^4) are suitably chosen for 
a given value of w. For instance, for w = —\ (the cosmo- 
logical constant) it is : 

{h,l2,h,h) = (2.6693,0.5950,0.0719,-3.0099) . 

At this point, one can wonder if it is possible to improve 
such a result considering even the radiation, although 
energetically negligible. Rather than inserting radiation 
in the H51|l . a more general approach in this sense is to 
consider the Hubble parameter descending from a unified 
model like those discussed in [ig. In such a scheme one 
takes into account energy density which scales as : 



p{z) = A [1 + 



having defined : 



1 



l + Zs 



P-a 



1 



Zb 



Zs = 1/s - 1 



Zb 



1/fe-l 



(57) 



(58) 



This model, with the choice (a,/3) = (3,4), is able to 
mimic a universe undergoing first a radiation dominated 
era (for z 3> Zg), then a matter dominated phase (for 
Zb -^ z <^ Zg) and finally approaching a de Sitter phase 
with constant energy. In other words it works in the way 
we are asking for. In such a case, the Hubble parame- 
ter can be written, in natural units, as H — y 3^ and 
one can perform the same calculation as in the quiessence 
case. As final result, it is again possible to find out a suit- 
able /(-R)-gravity model which, for numerical reasons, it 
is preferable to interpolate as f{R)/R: 



fJR) 

R 



1.02 X 



-0.69 X ( 



R 

Rt) 



R 

Rn 



-0.53 



-0.04 X (A)0.3i 

Rq 



X In(^) 



(59) 



where Rq is a normalization constant. This result once 
more confutes issues addressing modified gravity as in- 
compatible with structure formation prescriptions. In 
fact, also in this case, it is straightforward to show that 
a phase of ordinary matter (radiation and dust) domina- 
tion can be obtained and it is followed by an accelerated 
phase. 

4. In this paper we have discussed the possibility, in 
the framework of /(i?)-gravity models, to obtain a stan- 
dard matter phase followed by accelerated expansion. We 
have faced this task following two different lines: in a 
first case, assuming a gravity Lagrangian oc R ~ cii?~" 
we have outlined the mapping of f{R) gravity from the 
Jordan Frame to the Einstein one by means of a con- 
formal transformation. We have showed that, as mat- 
ter of fact, dust matter domination is possible in this 
approach and a late time accelerating phase naturally 
arises. When a minimally coupled scalar field Lagrangian 
is mapped back into the Jordan frame, one can deduce 
that modified gravity can imply an effective fluid contri- 
bution whose origin is the geometrical counterpart. The 
consequent EoS H31|) shows a constant barotropic index 
which, in relation to the choice of the model parame- 
ter n, can provide behaviours compatible with ordinary 
matter evolutions. Furthermore, such a parameterization 
allows also to describe phantom-like regimes, which turn 
out to be inspired by /(i?)-gravity. When a mapping 
into the Einstein frame is considered, this solution does 
not give iJ > so that the phantom - like behaviour is 
lost. This means that the conformal transformation does 
not necessarily imply physical equivalence of solutions 
and this difference among the two frames could repre- 
sent a pathology of back mapping conformal frame solu- 
tions. The debate on this point is recently growing and 
several studies have been pursued in this direction (see 
e.g. [33|)- In order to avoid pathologies related to the 
conformal transformation of solutions, another approach 
can be pursued, discussing the viability of dust matter 



regimes directly in the Jordan Frame. We have exploited 
a reconstruction procedure of /(i?)-gravity by means of 
the phenomenological (or observational) knowledge of the 
Hubble parameter in term of the redshift. The H(z) rate 
can be inspired by data fitting of speculative models like 
Q-essence or by means of unified schemes considering an 
ordinary matter-dominated expansion followed by a late 
time acceleration. Since in both cases /(i?)-gravity is 
compatible with the assumed H{z), one can conclude 
that modified gravity is sufficiently versatile to furnish 
the dynamics of a viable cosmic evolution. 

In summary, it seems that /(i?)-models are viable and 
cannot be ruled out only assuming mathematical equiv- 
alence of conformal transformations. However, consider- 
ing simple f{R) = foR"' models cannot completely solve 
the cosmic evolution problem since matter-dominated so- 
lutions for n y^ 1 are generally stable and do not give 
rise to late-time transitions to accelerated behaviors (see 
also ^3I)- The situation is similar to that faced several 
times in inflationary cosmology where graceful exit from 
infiation regime can result problematic for simple mod- 
els. In other words, f{R) gravity could be the way to 
bypass shortcomings as dark energy and dark matter at 
cosmological and astrophysical scales taking into account 
only o&serued ingredients as baryons, radiation and grav- 
ity but more realistic Lagrangians have to be taken into 
account than simple power-law ones. 

In the following appendix, we show that it is always 
possible to work out a /(i?)-model capable of match- 
ing with cosmological radiation, matter and accelerating 
epochs in succession. 

Appendix. Let us consider the reconstruction of mod- 
ified gravity. That is, first we consider the proper Hub- 
ble rate H , which describes the evolution of the universe, 
with radiation dominance, matter dominance, and accel- 
erating expansion. It turns out that one can find f{R)- 
theory realizing such a cosmology (with or without mat- 
ter) . The construction is not explicit and we need to solve 
the second order differential equation and algebraic equa- 
tion. It shows, however, that, at least, in principle, we 
could obtain any cosmology by properly reconstructing a 
function f{R) on theoretical level. 

The starting action is 



S^ 



d^x{P(0)i?H-g((/.)+Anatter} 



(60) 



Here P and Q are proper functions of the scalar field (j) 
and /Cmattor is the matter Lagrangian density. Since the 
scalar field does not have a kinetic term, we may regard 
(j) as an axilliary field. In fact, by the variation oi (f>, it 
follows 



O-P'(0)i? + Q'(0) 
which may be solved with respect to ( 



(61) 



(62) 



By substituting H62() into H60|l . one obtains /(i?)-gravity: 

S = Jd^x{f{R) + C,^,ttoA , 
/(i?) = Pi^{R))R + Q{^iR)) . (63) 

By the variation of the action (|60() with respect to the 
metric g^^, we obtain 

= -lg^APi4')R + Qm-Rf.^PW 

+V^V,P(</>) - 5/..V2p(0) + ip^, . (64) 
FRW equations are 

= -6g^P(0) - Q{c^) - QH '^^^f}^^^ + p , (65) 



dt 



= (4P- + 6ff2) P(0) + Q(0) 

d^pm)) , ,^dP{m) 



+2- 



dt 



AH- 



dt 



+ p. (66) 



By combining H64|) and H65|) and deleting Q {(/)), we find 
the following equation 



= 2 



d'Pm)) ^dP{4>{t)) 



dt^ 



2H- 



'AHP{(f>)+p + p. (67) 



As one can redefine the scalar field properly, we may 
choose 



= t 



(68) 



We also assume, as matter, for example, a combination 
of the radiation and dust: 

p = proa + pdoa , p ^ —a . (69) 

Here Pro and pdo are constants. If the scale factor a is 
given by a proper function g(t) as 



a = aoC' 



ait) 



(70) 



with a constant ao, Ea. (|66|l reduces to the second rank 
differential equation: 



= 2 



d^Pjct)) 
d(f>^ 



- 4g"(0)P(0) - 25'(0) 



u..dPi<f')) 



dt 



+ ^P,.o«o*e-4«(*)+p,oao^e-3«(^'. (71) 

In principle, by solving H71|l we find the form of P{(l))- 
Then by using (|65|l (or equivalently (|66|l ). we also find 
the form of Q{(f>) as 

Q{<P) = -6(5'(0))'P(0)-65'(0)^^ 

+p,oao ^e-49('A) + p,oao-3e-39(« . (72) 



Then in principle, any cosmology expressed as (|70|) can 
be realized by /(i?)-gravity. 

For example, we consider the case 



51 



9'{4>) = 9o + — , 



(73) 



Here Pq is a constant. Then by using H61|l . we find 



6.91 (.91 + g) (g - 2) 
aR 



(80) 



and without matter p 
Ea. (|71|l reduces as 



for simplicity. Then ^^ich gives 



= 



d^P 



■90 + -r 



gi\ dP 2(71 



P 



(74) 



whose solutions are given by the Kummer function (hy- 
pergeometric function of confluent type) as 

P = z''FKia,r,z), z^-^FK{a--f+l,2~r,z). (75) 



Here 



z = 909 



_l + gi± y/g'f + 2^1 + 9 



7=1± ^ , 

and the Kummer function is defined by 



i^K(a,7;z) 



^ a{a + 1) • • • (a + n — 1) z" 

;^7(7 + i)---(7 + "-i)^ 



Eq.{7nj) tells that the Hubble rate H is given by 



^ = ^" + T 



(76) 



(77) 



(78) 



Then when t is small, as H ^ 91 /t, the universe is 
filled with a perfect fluid with the EoS parameter w = 
— 1 + 2/3gi. On the other hand when t is large H ap- 
proaches to constant H — j go ^nd the universe looks as 
deSitter space. This shows the possibility of the tran- 
sition from matter dominated phase to the accelerating 
phase. Similarly, one can construct modified gravity ac- 
tion describing other epochs bearing in mind that form of 
the modified gravity action is different at different epochs 
(for instance, in inflationary epoch it is different from the 
form at late-time universe). 

We now investigate the asymptotic forms of f{R) in 
(|63|l corresponding to (|73|l . When (p and therefore t are 
small, we find 



P ^ Por , Q - -6P0.91 (.91 + «) 0' 



a-2 



(79) 



/(i?) ^ -^ (^^i(^i±i^ii^^^r^-f . (81) 



a — 2 [ a 

On the other hand, when (j) and therefore t are positive 
and large, one gets 



P 



Po'/'^""'^e9«'^ I 1 + 



(l-a)(7-a) 



-9 + 12a - 57 - 2a7 + 2a'^ 



X 1 + 



29o(t> 



-l + 9a-l^ 



90 






(82) 



Here Pn is a constant. Then we find 



/(i?) ^ l2glP„ 



.90 



Vl2.9o' 



9^7 

-2+9a--7 

-2a+7+l 



2a— 7 



X exp 



-|+9a-|7' 
-R^^l 



(83) 



This shows the principal possibility of unification of 
matter-dominated phase, transition to acceleration and 
late-time speed up. 
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